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^ ■ Abstract 

O ; 

^ ! We present an approach of computing the intersection curve C of two rational paramet- 
j^ I ric surface Si(m,s) and S2(w,t), one being projectable and hence can easily be imphcit- 
ized. Plugging the parametric surface to the implicit surface yields a plane algebraic curve 
G{v,t) = 0. By analyzing the topology graph ^ of G{v,t) = and the singular points on 
the intersection curve C we associate a space topology graph to C, which is homeomorphic 
to C and therefore leads us to an approximation for C in a given precision. 
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1. Introduction 



> 

(N 

'sj" , Computing the intersection curve of two surfaces is widely studied in CAGD [1, 16, 

19], which is popularly applied in CAD/CAM and manufacturing. Existing approaches 
can typically be classified into numerical and algebraic categories. A very recent method 



o 

^ ■ proposed by [10, 11], they implemented the typical process to get the numerical intersection 
curve. However, some important geometric features of the intersection curve, such as self- 
intersected points and cusps, might be lost due to their numerical approximation nature. 

S>^ I An efficient algorithm that is robust, accurate, and requiring the least user intervention is 
H I therefore needed. 

" ■ ■ A projectable surface can be regarded as a planar curve in an extended coefficient filed. 

We show that projectable surfaces can easily be implicitized by simple successive resultant 
computations. The projectable surfaces include many widely used surfaces, such as whirled 
surfaces, ruled surfaces and obit-based surfaces. The projectable surface is considerable 
based on the fact that there are significant advances on intersection computing for these 
modeling surfaces, such as quadrics [18, 20, 25], ringed surfaces [14] and ruled surfaces [10, 
13, 3]. 
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In this paper we compute the intersection curve of two rational parametric surfaces 
Si(m, s) and S2{v,t), with Si{u,s) being projectable whose imphcit form F{x,y,z) = 
is computed by our proposed method. By substituting the parametric surface S2(m,s) to 
F{x,y,z) = we get a plane algebraic curve G{v,t) = 0. By analyzing the topology graph 
'^ of G{v,t) = and the singular points on the intersection curve C we associate a space 
topology graph to C, which is homeomorphic to C and therefore leads us to an approximation 
for C in a given precision. 

According to the above process, one important work is the implicitization of a given 
surface. Many methods such as the resultant method [10, 21], the Groebner basis method [7] 
and the //-basis method [2] were proposed. The resultant method is comparable in the 
complexity of computation but not complete for general surfaces. The latter two methods 
are complete but not efficient in practical implements. It is still a problem to implicitize 
a general surface efficiently. For the projectable surfaces, we introduce an implicitization 
method using simple successive resultant computations. The method is more efficient than 
others in these special cases and it can be introduced to numerical computation. 

Another significant work is topology determination of a real algebraic curve. Existing 
methods of determining the topology are referred to [15, 12, 22]. Based on the concept of 
segregating box in [4] and real roots isolation of triangular system isolation [5], we propose a 
method to compute the topology graph of a planar algebraic curve inside a box. The given 
curve need not to be in a generic position. 

Some subtle discussions are proposed to refine the topology graph, since there are the 
points which make the correspondence between the plane curve and the intersection curve 
not one-to-one. These points are figured out and added to the topology graph. Then each 
edge of the refined topology graph is homeomorphic to its corresponding curve segment of 
the intersection curve. 

The rest of the paper is organized as follows. In Section 2, we give some notations. The 
implicitization method is also proposed for projectable surfaces. In Section 3 we outline 
the process to determine the topology graph of a plane curve. In Section 4, we refine the 
topology graph and compute the space topology graph of the intersection curve. In Section 
5, we give an algorithm to approximate the intersection curve. Some experiments are shown 
in Section 6 and we draw a conclusion in Section 7. 

2. Implicitization of rational projectable surfaces 

A rational parametric surface S{u,s) is defined by 

/ ^ e^ ^ fPii^^^) P2{u,s) P3{u,s)\ 

{x,y,z) = S{u,s)= I— -,— -,— , 2.1 

\qi{u,s) q2{u,s) q3[u,s)J 

where Pi,qi G Q[ti, s] and gcd{pi,qi) = l,i = 1,2,3. We assume that the parametrization 
is proper [17] and presents a non-degenerate surface. Since there still lacks an efficient 
implicitization method that applies to general parametric surfaces, we suppose one of the 



two parametric surfaces, say Si{u, s), to be a rational projectable surface that takes the form 

{x,y,z) = S{u,s)= {— -, — r,^rT • 2.2 

The surface (2.2) can be treated as a collection of the following planar curves with specified 
parameter s: 

(.,,) = S(.;.)=(^.^). (2.3) 

\qi{u]s) q2{u;s)J 

Since gcd{pi,qi) = 1 and max{deg(pi), deg(gi)} >l,i = 1,2, 

Res(gix - pi, q2y - P2, u) = l{s)L{x, y, s) (2.4) 

is not identically zero, where l{s) G Q[s] is the content of the resultant w.r.t. parameters 
x,y. Hence, L[x,y,s) is the primitive part of the resultant. 

If degj,(L) = 0, the surface (2.2) is a cylindrical surface over the a:|/-plane with the irre- 
ducible implicit equation L{x, y) = 0. To determine whether a rational surface is cylindrical 
over the coordinate plane is not hard, hence we consider only the non-degenerate case with 
deg^(L) > 1. We have the following theorem. 

Theorem 2.1. Let S{u, s) be a rational projective surface of the form (2.2) and its implicit 
equation is F{x,y,z) ^ Q[x,|/]. Then up to a constant multiple 

Fix, y, z) = Res{q3{s)z - ^3(5), L{x, y, s),s). (2.5) 

This theorem is a simplified version of theorem 2 in [23]. This implicitization method is more 
efficient than most existing approaches. Readers are referred to [23] for further details. 

2.1. Ruled surfaces 

A ruled surface is formed by a one-parameter family of straight lines moving along a 
curve, where the curve is the directrix, and the straight lines are called rulings. Ruled 
surfaces are widely used in geometric modeling, see [13, 10, 3] for related work. Precisely, a 
rational parametric ruled surface is given by 

C('^ g\ — f ao{u)+ai{u)s bo{u)+bi{u)s co{u)+ci{u)s \ /2 r-\ 

\ ^ ) y di(u) ' d2{u) ^ d-i{u) J \ • ) 

where ao;i, 6o;i, co;i, 7^ (ii;2;3 ^ Q[m]- We assume that the parametric equations in (2.6) are 
irreducible fractions, for instance, gcd(ao, ai,(io) = 1- Although the parametrization (2.6) 
does not agrees with (2.2), the following lemma allows us to transform (2.6) to a rational 
projectable surface. 

Lemma 2.2. Let S(u, s) he a rational ruled surface of the form (2.6). Then by a birational 
parameter transformation, S{u, s) can be reparameterized by 

g('^ g\ _ f ao(u)+ai{u)s bo{u}+bi(u)s - 
^ ' '^ y di{u) ' d2{u) ' 

where u, s are new parameters and ao-i,bo-i, Cq-i, 7^ Ji;2 G Q[u]. 
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Proof. Since S is a rational ruled surface, ai{u), bi{u), ci{u) can not simultaneously be 
identical to zero. Without loss of generality, we assume ci{u) ^ 0. By introducing 



we get 



_ 


co{u) + ci{u)s 






4(m) 


- u, 


s = 


4(m)s - co(n) 


= M, 



(2.7) 



which when substituted into (2.6) yields 
S(u, s) 

di{u) ' d2{u) ' dsiu) 

ao(u)ci{u)-ai(u)co{u)+ai(u)d3{u)s bo{u)ci{u)—bi(u)co{u)+b-i{u)d3{u)s - 
ci{u)di(u) ' ci(u)d2{u) ' 

ao{u)+ai{u)s bo(u}+bi(u)s - 
di(u) ' d2(M) ' 

This gives the parametrization of a projectable surface. D 

Once the ruled surface (2.6) is reparameterized by Lemma 2.2, we can apply Theorem 2.1 
to compute its implicit equation. 

Example 2.3. Let {x,y,z) = S{u,s) be a ruled surface given by 

( l-v?-2su 2u+s{l-v?) 

Since the parametrization is already projectable, we directly apply Theorem 2.1. First we 
compute 

L(x, y; s) = 4 1/2 + 4 s V + 4 s^y"^ - 8 s^ - 4 s^ + 4 x^ - 4. 

By removing the content (4^^ + 4), i.e., the gcd of the coefficients of L{x,y; s), we get the 
primitive part — s^ + y^ — 1 + x^. Then the implicit equation of the ruled surface is 

-z'^ +y^ -1 + x^ = 0. 

For comparison, readers can see [10] for the implicitization of the same ruled surface by 
computing the gcd of three resultants. 

2.2. Generalized revolution surfaces 

Revolution surfaces are also popularly used in manufacturing, such as porcelain modeling. 
A rational generalized revolution surface is defined by 



S{u,s) = (^ 



Pijs) 2u P2{s) l-u^ Pijs) 
qi{s) l+«2' g^is) 1+^2' q3{s) I ' 



When pi/qi = ^2/92? this defines a usual revolution surface rotating around the z-axis. 
Since this parametrization agrees with the form (2.2), we can use Theorem 2.1 directly for 
implicitization. 
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2.3. Orbit-based surfaces 

An orbit-based surface is a rational surface formed by translating a plane curve with 
its posture unchanged along a space curve. One can find that the orbit-based surface is a 
special case of the sweep surface. For instance, a tube surface can be defined by a circle set 
whose center follows a space curve 



) 



C/',, „^ _ { 2u _ pijs) l-u' _ P2{s) P3is) \ 
'^ya,i,) yi+u'^ gi(s)'l+«2 q2{s)^q3{s)J 

where {pi/qi,P2/l2^'Pz/l3) represents a space curve, see Example 6.2 in Section 6. 

Remark 2.4. We have shown an efficient approach of implicitizing rational pro jectable sur- 
faces (up to a birational parameter transformation) . Efficiently implicitizing arbitrary ratio- 
nal parametrized surfaces is still left an open problem. Notably the following process applies 
to compute the intersection loci of two general surfaces with one rational parametrized and 
the other being in implicit form. 

3. Topology determination of planar algebraic curves 

By the method proposed in Section 2, we compute the implicit equation F[x, y,z) = 
of the projectable surface Si. Substituting 82(1*, t) to F{x, y, z) = yields a plane algebraic 
curve G{v,t) = 0. We next determine the topology graph ^ of the curve G{v,t) = inside 
a given rectangle. 

There are many related work about computing the topology of algebraic curves [6, 8, 12, 
15]. We prefer the methods which need not require the curve to be in a generic position and 
need not to compute a Sturm-like polynomial sequence. We use the concept of segregating 
boxes in [4] to determine the adjacency relationship when we compute the topology of 
algebraic curves and real roots isolation of triangular system [5] to get the critical points of 
the curve. We will compute the topology of the curve inside a bounding box. 

Definition 3.1. A point Pq = (fo,to) is said to be a singular point on the curve G{v,t) = 
if G{vo,to) = G„(fo,to) = Gt{vo,tQ) = 0. ^4 point Pq = (fo,to) is said to be an v-critical 
point (resp. t-critical point) of C if GiyQ.to) = and Gt{vQ,tQ) = (resp. Gy{vQ,tQ) = 0). 

Definition 3.2. Let P be a point on the curve G{v,t) = 0. The left (right) branch number 
of P is the number of curve segments of C that passes through P from the left(right) in a 
small neighbor of P. 

The following definition is taken from [4]. 

Definition 3.3. Let f{x,y) G Q[x, y] be the defining polynomial of an algebraic curve and 
g{x) its discriminant with respect to y. A Segregating box of an x-critical point P : (a,/3) 
of an algebraic curve f{x,y) = such that g{a) = f{a,f3) = is a rectangle [a,b] x [c,d\ 
containing P inside such that 

1. There is no real roots of g{x) = in [a,b]. 
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2. There is no real roots of g{x) = f{x,y) = in [a,b] x [c,d]. 

3. The upper and bottom boundaries of [a, b] x [c, d] have no intersection with the curve 
fix,y) = 0. 

Suppose G{v,t) is square free and contains no univariate factor(s) in v^. 

We will compute the topology of curve inside a box B = [A, B] x [C, D] G Q^, and then 
determine the topology graph ^ = {V,S}, where V and S inside the box are defined as 
follows. 

• P is a set of points in the v — t plane: 

V = {Pi,j = («i, A,i), 0<t<s,0<j<Si} (3.1) 

where s,Si G N and (oj, f3ij) are towers of real algebraic numbers such that ao < «i < 
■ ■ ■ < as and /Sj^o < A,! < ■ ■ ■ < A^sr The points Pij shall later be solved from the 
triangular systems Sj = {hi{v),gi{v,t)} and then represented by the isolation boxes 
Bij = [ai, bi] X [cjj, dij]. Note that ^ or "P has s + 1 columns of points. 

• 8 = {(Pi, P-2)\Pu P2 e P, s.t. either Pi = P,,^, P^ = P,+i,, or Pi = P,,^, P^ = P.,p+i}. 
In the first case, the edge is called non-vertical, while in the second case, the edge is 
called vertical. We shall further assume no intersection between any two edges except 
at the endpoints. 

The following process outlines our approach to computing the topology graph ^: 

Algorithm 3.4. Compute the topology of a planar algebraic curve G{v,t) = inside a 
bounding box B = [A, B] x [C, D] . 

Step 1 Compute d{v) := {v - A){v - B)G{v,C)G{v,D)Res{G,^,t). 

Step 2 Solve for the real roots of the triangular system S = {d{v), G{v, t)} = by the real 
root isolation method given in [5]. We take only the real roots inside B. 

Step 3 For each critical point P = (aj, /3ij) solved from Step 2, construct a segregating box 
[aj,6j] X [cij,dij]. The number of the left branches of P is the number of roots of 
G{ai,t) = inside the interval [cij,dij]. Note that the line segments {cj < v < bi,t = 
Cij, dij} have no intersections with the given curve. The number of right branches for 
each critical points is similarly computed. 

Step 4 Construct a topology graph of G = inside B. 



^ The topology of a given curve C : G{v, i) = is the same as that of the curve defined by the square 
free part of G{v,t). Moreover, if C contains vertical line(s), these vertical lines can later be added after the 
main part of C is analyzed. 



4. Topology determination of the intersection curve C 

Geometric character points on the surfaces Si and S2, such as ordinary irregular points, 
cusps and self-intersected points, are crucial to determining the space topology graph home- 
omorphic to the space intersection curve C. 



Definition 4.1. Let S{v,t) be a rational parametric surface. A parameter (fo,to) is said to 

f """ f 



be regular if ^ and ^ are linearly independent at {vo,to); otherwise (fo,to) is said to be 



irregular. 

Definition 4.2. Let F{x,y,z) = be an implicit surface (of a parametric surface S{v,t)). 
A point Q = {xo,yo,zo) is a singular point on F{x,y,z) = if Xo,yo,Zo is a solution for 
F = ^ = ^ = ^ = 0; otherwise Q is said to be non-singular. 

Suppose a projectable surface Si (up to a birational parameter transformation) and a 
rational surface S2 are given by 



Si(m, s) 
S2{v,t) 



fl{u,s) f2iu,s) fsjs) 

gi{u,s)' g2{u,s)' g3{s)y .^ ^^ 

Pi{v,t) P2{v,t) paCM)^ ^ ■ ' 

qi{v,t)' q2{v,t)' qsiv^t) 



and F{x,y,z) = is the implicit equation of Si computed by Theorem 2.1^. If a point 
S2{v,t) lies on the intersection curve C, the parameter (f,t) must satisfy F{S2{v,t)) = 0. 
Let 

G{v,t) := F{pi{v,t)/qi{v,t),p2{v,t)/q2{v,t),ps{v,t)/q3{v,t)), 

and define the square-free part of the numerator by 

G{v,t) = Sqf(numer(G(t;,t)) = 0. 

The intersection curve C of the surfaces Si and S2 is then determined by 

G{v,t) = 
^7 . ; x=pi{v,t)/qi{v,t) .^2) 

' ^ y = P2{v,t)/q2{v,t) 
z =p3{v,t)/q3{v,t). 

Based on the topology of plane curve G{v,t) = 0, one can divide the intersection curve C^ 
into different curve segments. We now consider the parameter correspondence between the 
plane curve G{v,t) = and the intersection curve C^ of the two surfaces. 



^During implicitization, there might be some real pomts on F{x,y,z) ~ but not on the parametric 
surface Si(m, s). A connected component formed by such points is caUed a geometric extraneous compo- 
nent [9, 11]. Such points are singular points of the surface F{x, y, z) = 0. As mentioned in [11], we need to 
check singular points in order to verify if they belong to a geometric extraneous component, which must be 
removed. 



4.I. Self-intersection points 

Lemma 4.3. Except for the singular points on the plane curve G{v,t) = and the irregular 
parameters for the surface S2, the tangent of the point on the intersection curve (4.2) is 
different from zero. 

Proof. In the neighborhood of a point {v,t), one can regard v is the function of t, i.e., 
V = v{t) such that G{v(t),t) = 0, and S2{v,t) = S2{v{t),t). By Imphcit Function Theorem, 
one has ^ = ~^/^- Hence the tangent vector w to the intersection curve C^ at the point 
S2(w,t) is 

_dS2{v{t),t) _dS2dv{t) as2 

dt dv dt dt 

Since, then up to a constant multiple 

aSa dG 9S2 dG 
w = . 

dv dt dt dv 

For a parameter pair (vo,to), if neither ^ x ^ nor (^, ^) vanish at (uo,to)) i-e-) the 
parameter (wo,to) is regular on the surface S2 and the point (fo,to) is non-singular on the 
plane curve G{v, t) = 0, we have w 7^ 0, which means that the tangent of the point S2(t'o, ^0) 
is different from zero on the intersection curve C^ . D 

According to Lemma 4.3, besides those singular points {v, t) on the plane curve G{v, t) = 
(which are already computed in Section 3), the irregular points on the surface S2(f , t) shall 
also lead to singular points on the intersection curve C^ . These irregular points {v,t) are 
solutions for 

Giv,t) = 0. ^ '^^ 

A bad situation may occur for (4.3), that is, the surfaces S2 has a irregular parameter locus 
that corresponds to points on the intersection curve C^ . Under this situation, (4.3) has an 
infinite number of solutions. We shall first remove this common irregular parameter locus 
which shall be treated as a special curve component. 

As known that a point Q is a cusp of the C^ if the tangent of Q is vanish. Hence, the 
parameters corresponding to the cusp are included in the singular points of G{v,t) = and 
the solutions of (4.3). 

Notably, some self-intersection points of the intersection curve C^ may neither correspond 
to singular points on the plane curve G{v, t) = nor correspond to the irregular parameters 
on the surface S2{v,t). See the following example. 



Example 


4.4. 


Given two 


surfaces 






Si: 


H: 


-- {su^ + s 

■ su 

■ u 



l)/2 




l)t-l/Av 
l){t-l)t 




Figure 1: Intersection curve of Si and S2 

where Si is an elliptic paraboloid with whose implicit equation is 2x — y^ — z^ — z = 0. 
By computation G{v,t) = 2v + t'^ + t^ + t^v/2 - 2t^ - tv/2 - v^/lQ -t'^ + t = Q, which 
has no singular points. However, the intersection curve C^ has two self-intersection points 
(0,0,0) and (32,-8,0) whose corresponding parameters on G{y,t) = are {(0, 0), (0, 1)} 
and {(32, 0), (32, 1)} respectively (Figure 1). One can check that these two points are regular 
in S2(f ,t). 

This subtle situation did not draw enough attention in previous work, such as in [10]. 
This might omit some self- intersection points on the curve C^, whose corresponding param- 
eters are both non-singular points on G{v,t) = and regular points on surface S2(f,t). 
See Figure 2. Therefore, to ensure that the associated space graph ^^ mapped from ^ is 
homeomorphic to the intersection curve C^ , here we refine the topology graph ^ by adding 
the following character points of the intersection curve C^: 

{Pij = {ai, Pij) I {x{ai, Pij), y{ai, Pij), z{ai, Aj)) is a self-intersection point on C^} (4.4) 

with < i < s,0 < j < Si- Note that this might reintroduce some points {v,t) that are 
already computed in Section 3. Our principal here is not to omit any possible {v,t) that 
is crucial to the topology of the space curve C but do not pursue the complement of the 
previous computed {v,t). 




/\ 


/ 


1 ^ ^ 

' 1 , 






/ 



Figure 2: Self-intersection point overlooked by numerical approximation 

We now show how to compute the self- intersection points of (4.2). The self-intersection 
points {x,y,z) of C^ are those self-intersection points on S2(f,t) that also lie on Si, and 
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hence can be solved through: 

( x{v,t) — x{u, s) = 

I y{v,t)-y{u,s) = , . 

z{v,t)-z{u,s)=0 ^ ' 

G{v,t) = 0,G(m,s) = 0. 

Removing the denominators yields 

)qi{u, s) - pi{u, s)qi{v, t) = 

(4.6) 



' pi{v,t)qi{u,s) -pi{u,s)qi{v,t) = 
P2{v, t)q2{u, s) - p2{u, s)q2{v, t) = 
Psiv, t)q3{u, s) - p-iiu, s)q-i{v, t) = 

^ G{v,t) =0,G(m,s) = 0. 



Except for the trivial solution set {v = u^t = s}, the system (4.6) can only has a finite num- 
ber of solutions which correspond to the self-intersection points of the intersection curve 
C^ . These solutions can be obtained by Ritt-Wu's characteristic set method [26], supported 
by the maple packages Wslove and Charsets by D.K.Wang and D.M. Wang respectively. 
The packages can be downloaded from http://www.mmrc.iss.ac.cn/~dwang/soft.htm 
and http : //www-calf or . Iip6 . f r/~ waing/epsilon/. 

Notably during the zero decomposition in Ritt-Wu's method, some multiplicities of the 
solutions may be lost. However, we do not care about these multiplicities, since the multiple 
solutions are corresponding to the cusp points, which are already computed. 

Theorem 4.5. Except for the singular points on the plane curve G{v,t) = 0, the irregu- 
lar parameters of the surface S2 from (4.3) and the parameters of self-intersection points 
from (4.6), there is a one-to-one correspondence between the plane points on the plane curve 
G{v,t) = and the space points on the intersection curve C^ of the two surfaces. 

Proof. By Lemma 4.3, the tangent is well defined at the space point on the intersection 
curve C^ of the two surfaces. Furthermore, the self-intersection points are excluded in 
equation system (4.6). Hence, except for these points, there is a one to one correspondence 
between the plane curve and the intersection curve of the two surfaces. D 

4-2. Space topology graph 

We now determine the space topology graph W^ , whose vertices are mapped from the 
vertices of the refined graph of W: 

{P-j = ((x(ai, (3ij),y{ai, (3ij), z{ai, Aj)) : a^, /3jj)}, < i < s,0 < j < Si, 

and edges (-Pj^jji -Pjaj^) I'epresent the topological connections between the points -Pj^^j and 
^ti 72 °^ ^^ ■ ^o^^ ^J^at any two edges {Pl^P^) and {P^^Pl) can only share one of their 
endpoints. 

Algorithm 4.6. Compute the refined topology graph ^ of the irreducible curve G{v,t) = 
and the space topology graph ^^ of the space curve (4.2). 
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1. Compute the critical and singular points of G{v,t) = and determine the topology 
graph W by Section 3. 

2. Compute the cusp and self-intersected points of the space curve (4.2) by the methods 
in Section 4-1- 

3. Refine the the topology graph ^ by adding the parameters corresponding to the points 
computed in step 2, as well as the edges. 

4. Map the vertices of the refined graph of ^ to space points to get the vertices of the 
space graph '^^ , and correspondingly map the edges in the refined graph of ^ to the 
edges of the space graph ^^. 

5. // any pair of the edges [Pi , Pj) and (P/, P/) in ^^ have no intersection points except 
at the endpoints, then output '^^ and'^^ . Otherwise, if two edges {Pl,Pj) and {Pl,Pl) 
intersects at transver sally, add subdivision vertices P*^ = (i;*^,t^^), m = 1, . . ., niQ to 
$f between P/ and Pj , P2^ = {v2^,t2^),n = 1, . . . , kq between P/ and P^ , such that the 
line segments {Pi Pl')XPu,Pui--; (PiJ^P/) and {P^P^^, {P*/ ,P*/), . . ., 
(P^ , Pj ) are topology edges ofi^^. Co to Step 4- 

In Step 4 and 5, one can determine the intersection of two line segments (P/, Pj) and 
{Pi, Pi) using the bracket formulas in [3]. 

Lemma 4.7. In Algorithm 4-6, the output space topology graph ^^ is homeomorphic to the 
intersection curve C^ . 

Proof. For two points Pj, Pj on graph ^, their edge (Pj, Pj) corresponds to a curve segment 
PiPj of G{v, t) = 0. Map them to the space topology graph ^^ , we get P/, P^ and their line 

segment {Pl,Pj) on '^^ . The space curve segment Pi Pj on C^ is then equivalent to PiPj- 

If there exists no singular points on PiPj, then by Theorem 4.5, PlPj is a continuous curve 

segment. According to Algorithm 4.6, there has no cusp or self-intersection points on Pi Pj 

except for the endpoints. Hence PlPj is homeomorphic to the line segment (Pj , Pj ). 

If the edge (Pj , Pj ) intersects with another edge {Pk , Pi ), one can subdivide these two 
edges by adding a finite number (mo and uq) of points in Step 5. This is based on the facts 

that a curve segment can be approximated in any precision by line segments, and that PlPj 

and PlPl have no intersection points excepted in endpoints. Since any two line segments 
has no intersection point except in the endpoints, the line segments are then the edges of the 
space topology graph ^^ . Then the output ^^ is homeomorphic to the intersection curve 
C^. D 

Remark 4.8. In Algorithm 4-6 and Lemma 4-7, we simplify the discussion by assuming 
G{v,t) irreducible. In fact, the algorithm and lemma can be enhanced for general cases 
G{v,t) by factorizing G{v,t) to irreducible factors Gi{v,t) and decomposing the curve to 
the components in assumed form. However, to combine these decomposed components, we 
should compute the common points of Gi{v,t) = and add them to each topology graph ^i, 
since these intersections may be lost in the numerical computation. 

11 



Example 4.9. Continue with Example 4-4> the topology graph ^ of G{v,t) = not in- 
cluding the self-intersected points of C^ , its mapped topology graph ^^ and the numerical 
intersection are shown in Figure 3. '^^ and the numerical curve loss the self-intersected 
points of C^ . 




Figure 3: ^, ^^, numerical intersection and enlarged neighbor 

Adding two self-intersected points, we get the refined'^ , ^^ and the numerical intersection 
as Figure 4- The red diamond points are added since they correspond to the self-intersected 
points. The refined '^^ and the numerical curve have same topology with C^ . 




Figure 4: Refined ^, ^^, numerical intersection and enlarged neighbor 

The self-intersected points are preserved as a topology vertex in our numerical intersection 
now (right one of Figure 4) while it may be lost in a numerical approximation (right one of 
Figure 3). 

5. Approximation of the intersection curve 

Since the plane topology graph ^ and the space topology graph W^ are both determined 
in the last section, we now approximate the intersection curve C^ within a given precision e. 
Our principal is to subdivide the plane topology graph W, and hence the space topology 
graph '^^ is simultaneously subdivided. 

Suppose the vertices of the plane graph ^ are 

"^ = {Pijivi, U,j) e [ci, hi] X [cjj, (iij], < « < s, < j < Si}. 

Now we consider the boxes Bj = [fj,fj+i] x [C,D]. Suppose there are rrii edges of ^ 
enclosed in Bj, i.e, mi curve branches originate from the vertices {-Pi,j, < j < Sj} and end 
at {Pj_|.ij,0 < j < Sj+i}. Now for each i, rename the vertices of ^ on the vertical lines 
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V = Vi and v = Vi+i as {Lj,0 < j < m.j — 1} and {Rj,0 < j < rrii — I}, respectively. 
Note that some x-critical points may repeat several times in L (or R) list. The following 
procedure tells how to approximate the intersection curve C^ within a given precision e: 

Algorithm 5.1. Approximate C^ within a given precision e from the topology ofG{v, t) = 0. 

1. We first deal with the vertical edges in the graph ^ (note that the v— coordinates 
of these vertical edges are zeros of the content V{v) of G{v,t)) if they exist. For 
a vertical edge v = Vi in ^, its corresponding component in the space curve C^ is 
{x{vi,t),y{vi,t),z{vi,t)) = S2{vi,t). 

2. For i = 0, . . . , s — 1, execute the following steps. 

(a) Let L = Vi+i — Vi, and let A^, be the minimal integer larger than L/e; if Ni = 1, 
we set A'j = 2. Let ti^k be the roots set of G{vi + k^,t) = inside [C,D], 
I < k < Ni — 1. Arrange ti^k from bottom to up, we have ti^k = {^ik^ • • • > ^Tk~^}- 
Note that points Qlf^ := (f, + k^,tlf^) are not x-critical points of C. 

(b) Get the two lists {Lj, < j < nii — 1} , {Rj, < j < m,j — 1} as mentioned before. 

(c) For each planar point Lj,Rj or Q^f^ (sketch shown in Figure 5), denoted by 

{va, ta,i3), compute its corresponding spatial point Lj, i?j or Q^lj., i.e., 82(^0, ta^is)- 

(d) For each spatial line segment in : 

{{Lj,Q i,l),(Q i,fc,<5 i,k+l)l<k<N,-2,{Q i,N,-l^Rj)} 

{j = 0, . . . ,mi — 1) (assuming the endpoints of the line segment are Pi{xi,yi, zi) 
and P2{x2,y2, Z2)), check the the Hausdorff distance between the line segment 

[Pi, Pi) and its corresponding curve segment PlPl in C^ . If 

Dis{{Pl,Pl),PlPJ)<e, (5.1) 

does not hold, subdivide the planar graph ^ until all the mapped subdivided 
spatial graph ^^ segments satisfy the above condition. 

(e) For any pair of line segments (P/, Pj) and (P/, P/), we subdivide them to topol- 
ogy graph edges if they have intersection points. It means that any pair of edges 
(P/, pj) and (P/, P/) has no intersection points except for endpoints in the sub- 
divided topology graph. 

Remark 5.2. In Step 1, since we use isolation interval for Vi G [aj,6j],fj+i G [aj+i,6i+i], 
we let L = (flj+i -|- 6j+i — a.i — bi)/2 in practical computation. If [ai + hi)/2 + L/Ni < hi (or 
(Oi+i + 6i+i)/2 - L/Ni > di), refine [a^, k] (or [qj+i, k+i]). 

Consider the planar curve segment subdivided in t^t-plane, all the endpoints of the seg- 
ments are the topology vertices on the curve G{v,t) = 0, and the endpoints of the corre- 
sponding spatial curve segments are the topology vertices on the intersection curve. It is 
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Figure 5: The arrangement of the points in B^ = [vi, Wi+i] x [C, D] 



clear that the approximation precision between the intersection curve and the numerical 
approximation is less than e since all the curve segments satisfy (5.1). 

The exact Hausdorff distance in (5.1) is not easy to compute, we actually compute the 

numerical distance instead. Choosing m points from PfP-j, we then compute the max 
distance between these points and {Pl.P^) as the distance of (5.1). 

Theorem 5.3. The numerical intersection is homeomorphic to the intersection curve C^ 
and convergence to C^ in parameters subdivision process based on ^ . 

Proof. Since the numerical computation is a parameter subdivision process based on the 
refined topology graph ^, by step 2(e), it is actually a subdivision of the space topology 
graph ^^ . According to Lemma 4.7, the numerical intersection is homeomorphic to C^ . 

Since all the character points are computed as vertices, they are preserved in the param- 
eter subdivision process. We can find that the numerical intersection curve converges to C^ 
as the line segments approximation. D 

The above numerical intersection is a line segment approximation. For further consid- 
eration, we can give the B-spline approximate intersection. A method is proposed in [24] 
to approximate a give space curve based on its topology graph ^^. For each vertex of '^^ , 
we can compute the left and right tangent directions as well as the osculating planes. On 
an ordinary point, the left and right tangents and normal direction are consistent as well 

as the osculating planes. Consider a space curve segment PlP^ with the tangent direction 
and osculating planes at endpoints, we can construct the cubic Bezier curves to approximate 

P(P2- Then rewrite the Bezier spline curve to B-spline curve with proper knots selection. 

Comparing with the line segment approximation, there are at least three advantages 
in cubic B-spline approximation. The first one is that the approximation B-spline is C^ 
continues except in the cusps while the line segment approximation is only C° continues. 
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Then second is that the cusps are preserved. Finally, the number of approximate curves 
segments is much less than that in line approximation. 

6. Experiments 

To illustrate our algorithm, we will give some examples in this section. Some of them 
are taken form [13] and [10] for comparison. 

Example 6.1. Consider the intersection of a cone and an elliptic cylinder [13, 10] 



Si 

S2 



l-v?+s{l-v?) 2 u+2 si 
1+^2 ; 1+^2 

l-v"^ 2v+t{lW) -, I 
l+«;2 1 l+t,2 ) ^ "f 



l + S 



t 



Since Si is projectable, one can compute its implicit equation and get the {v,t)-plane curve 
equation 

G{v,t) =t{l+v^){v-l) = 0. 

The solution of line t = corresponds to the red circle in Fig. 6, which is the directrix of both 
surfaces. The line v = 1 corresponds to the common ruling (blue line in Fig. 6) where the 
two surfaces meet tangentially. The intersection point of the circle and the line is (0, 1, 1) 
which corresponds to the singular point (1,0) ofG{v,t) = 0. 




Figure 6: A Cone and an Elliptic Cylinder / Topology of the plane curve 

This example was involved in a more careful discussion in [13] and the results are divided 
to four parts. Comparing with [10], we add the intersection point (1,0) in the space '^^ . 
Otherwise, the numerical intersection may consist of two components separately: a numerical 
circle no passing through (0, 1, 1) and a numerical line passing through (0, 1, 1). 

Example 6.2. As mentioned in Section 2, we construct a tube surface 



Si = (s3 + s, 



s^+s'u'^+u l-u'^+2s+2su'^ 



l+«2 ! 2(l+n2) 



formed by a circle following the space curve {s^ + s, s^, s). And the surface S2 is a whirled 
surface without much restriction. 



Uv v{'i+v'^) 



C _ / _ 2(-l+f2)t) 
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The implicit equation of Si is 

F{x,y,z) = 

-1 + 256 yz'^ + 16y + 384 x^z^y - 768 xzy^ - 52 y^ 
+832 z'^x'^ + 208 z^ + 192 yS^ + 192 xzy + 64 x^ 
+64 ^6 + 512 y^z^ - 768 y^xz - 128 i/^x^ - 52 z^ 
-128 2/^2 _ 384 x^z + 256 y^ - 96 yx^ + 320 y^x^ 
- 128 2/3 - 768 yxz^ -16x2 + 64 x^ + 64 y^ 
+192 2/2^^ _ 768 z^x + 208 y'^ + 416 1/^^^ = 0. 

We omit the intersection equation G{v,t) = for its long expression. The topology of the 
plane curve and the numerical intersection curves (red curves) are illustrated as the following 
figures (See Fig. 7). 






Figure 7: Topology of the plane curve / Numerical intersection curves 



Example 6.3. The given surfaces are 




and S2 



Si is an elliptic paraboloid whose implicit equation is 

F{x, y, z) = y'^ + z'^ — 2x + z 



2, vt 



1+v^ 







and S2 is a cone. Then G{v,t) = Gi{v,t)G2{v,t) consists of two irreducible factors as 
Gi{v, t) = t and G2{v, t) = tv^ + 3 1'"' + 6 tv"^ + t + 2 vH — 1. It means that the intersection 
of Si and S2 consists of two components: 



Ci:{ 



' Gi(t;,t) = 

l+?)2 

2, vt 



X 

y 

z ■■ 



and C. 



l+v^ 



t 



G2iv,t)=0 

1+v^ 

2 "'^ 



l+v^ 



t 



C{ is an isolate point as (0, 0, 0) and C^ is a quartic space curve with a self-intersected point 
(0, 0, 0) with the parameters ±(\/3/3, 0). We can find that the point (0, 0, 0) plays different 
roles: 1) an isolate point of C{ corresponding to (t',0) at S2; 2) a self-intersected point of 
Cl; 3) the intersect point of €[ and C|; 4) the singular point of S2. 
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Figure 8: Numerical intersection curves 



7. Conclusions 



The paper proposes an improved method for approximating the intersection curve of two 
rational surfaces, one being projectable. For two given rational surfaces, we enhance some 
key steps in the typical process of intersection analysis. Our method is simpler in impliciti- 
zation and adapt to more surfaces, similar enhancement in planar topology determination. 
As another important improvement, we refine the topology graph by adding more singular 
points as well as their corresponding points of the intersection curve. Then the space topol- 
ogy graph is homeomorphic to the intersection curve. And the numerical curve based on the 
space topology graph converges to the intersection curve in subdivision process. Based on 
the refined topology graph, we can approximate the intersection curve with B-spline curve 
and we will give the more details in the further paper. 

Acknowledgement 

This work is partially supported by National Natural Science Foundation of China under 
Grant 10901163, 11001258, a National Key Basic Research Project of China (2011CB302400) 
and a President Fund of GUCAS. The authors also wish to thank the anonymous reviewers 
for their helpful comments and suggestions. 

References 

[1] Barnhill, Farin, G., Jordan, M. and Piper, B.R.: Surface/surface intersection. Computer Aided Geo- 
metric Deisgn, 4, 3-16, 1987. 

[2] Chen, F., Zheng. J. and Sederberg, T.W.: The mu-basis of a rational ruled surface. Computer Aided 
Geometric Design, 18, 61-72, 2001. 

[3] Chen Y., Shen L.Y. and Yuan, CM.: Collision and intersection detection of two ruled surfaces using 
bracket method . Computer Aided Geometric Design, 28, 114-126, 2011. 

[4] Cheng, J.S., Gao, X.S. and Li, J.: Ambient isotopic meshing of implicit algebraic surface with singu- 
larities, CASC'2009, LNGS 5743, Springer, 89-93, 2009. 

[5] Cheng, J.S., Gao, X.S. and Yap C.K.: Complete numerical isolation of real roots in zero-dimensional 
triangular systems. Journal of Symbolic Computation, 44, 768-785, 2009. 

[6] Cheng, J.S., Lazard, S., Penaranda, L., Pouget, M., Rouillier,F., and Tsigaridas, E.: On the topology of 
planar algebraic curves. In SCG '09: Proc. of the 25th Annual Symposium on Computational Geometry, 
pages 361-370, New York, NY, USA, 2009. ACM. 

[7] Cox, D.A., Little, J. and O'Shea, D.: Ideals, Varieites, and Algorithms. Springer, Berlin, 2nd ed., 1996. 

17 



[8] EmeliyanenkOjP., Kobel,A., Berberich, E. and Sagraloff, M.: Arrangement Computation for Planar 

Algebraic Curves. In Symbolic Numeric Com,putation (SNC), 2011. 
[9] Fioravanti, M., Gonzalez- Vega: On the geometrie extraneous components appearing when using im- 

plicitization. Mathematical Methods for Curves and Surfaces, Nashboro Press, 157-168, 2005. 
[10] Fioravanti, M., Gonzalez- Vega, L. and Necula, I.: Computing the intersction of two ruled surfaces by 

using a new algebraic approach. Journal of Symbolic Computation, 41, 1187-1205, 2006. 
[11] Fioravanti, M., Gonzalez- Vega, L. and Necula, I.: On the intersection with revolution and canal surfaces. 

Mathematics and Visualization: Algebraic Geometry and Geometric Modeling, 169-183, 2006. 
[12] Gonzalez- Vega, L., Necula, I.: Efficient topology determination of implicitly defined algebraic plane 

curves. Computer Aided Geometric Design, 19, 719-743, 2002. 
[13] Heo, H.S., Kim, M.S. and Fiber, G.: The intersction of two ruled surfaces. Computer- Aided Design, 

31, 33-55, 1999. 
[14] Heo, H.S., Hong, S.J., Kim, M.S., Fiber, G.: The intersection of two ringed surfaces. The Eighth Pacific 

Conference on Computer Graphics and Applications (PG'OO), 146-153, 2000. 
[15] Hong, H.: An efficient method for analyzing the topology of plane real algebraic curves. Mathematics 

and Computers in Simulation, 42, 571-582, 1996. 
[16] Hoschck, J., Lasser, D.: Fundamentals of Computer Aided Geometric Design. Wellcsley, MA: A.K. 

Peters, 1993. 
[17] Li, J. Shen, L.Y. and Gao, X.S.: Proper reparametrization of rational ruled surface. Journal of Com- 
puter Science and Technology, 23(2): 290-297, Mar. 2008. 
[18] Miller, J.R., Goldman, R.N.: Geometric algorithms for detecting and calculating all conic sections in 

the intersection of any two natural quadric surfaces. Graphical Models and Image Processing, 57, 55-66, 

1995. 
[19] Patrikalakis, N.M., Maekawa, T.: Intersection problems. In: Farin, G., Hoschek, J., Kim, M.S. (Eds.). 

Handbook of Computer Aided Geometric Design. North-Holland, Amsterdam, 623-649, 2002. 
[20] Sarraga, R.F.: Algebraic methods for intersections of quadric surfaces in GMSOLID. Computer Vision, 

Graphics and Image Processing, 2, 222-238, 1983. 
[21] Sederberg, T.W., Zheng, J.: Algebraic methods for computer aided geometric design. In: Farin, G., 

Hoschek, J., Kim, M.S. (Eds.). Handbook of Computer Aided Geometric Design. North-Holland, Ams- 
terdam, 363-387, 2002. 
[22] Seidel, R., Wolpert, N.: On the exact computation of the topology of real algebraic curves. SCG'05, 

Pisa, Italy. 107-115, 2005. 
[23] Shen, L.Y. and Yuan CM.: Imlicitzation using univariate resultants. Joural of System Science and 

Complexity, 23, 804-814, 2010. 
[24] Shen, L.Y., Yuan CM. and Gao X.S.: Certified approximation of parametric space curves with cubic 

B-spline curves. Submitted. 
[25] Wang. W.: Modelling and processing with quadric surfaces. In: Farin, C, Hoschek, J., Kim, M.S. 

(Eds.). Handbook of Computer Aided Geometric Design. North-Holland, Amsterdam, 777-795, 2002. 
[26] Wu, W.T. Basic Principle of Mechanical Theorem Proving in Geometries (in Chinese). Beijing: Science 

Press, 1984; English Version, Wien: Springer, 1994. 



18 



